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Computing parton distribution functions in QCD is a formidable task. It is in- 
trinsically non-perturbative and and thus hard to calculate from first principles. 
On the other hand, knowledge of these distribution functions is required in any 
dynamical process involving hadrons. This note describes a new approach to the 
problem. It is shown how this now approach leads to saturation of the gluon 
density at small x. 



1 Introduction 

In perturbation theory, one considers evolution of soft gluon density in x, or 
both. In the double log approximation of DGLAP formalism, terms of the form 
aln(a;) ln((5^) are summed. This is done by considering the so-called ladder 
diagrams where both longitudenal and transverse momenta are ordered. 

For very small values of x it may be more appropirate to sum terms of 
the form Q!ln(a:). This is done in BFKL formalism where one has ordering in 
longitudenal momenta while assuming that everything is happening at more or 
less the same transverse area. It is important to realize that both DGLAP and 
BFKL are evolution equations and as such require initial parton distributions 
as input. 

Regardless of whether BFKL or DGLAP is more appropirate for descrip- 
tion of soft gluons, they both have one outstanding feature in common; they 
both predict a sharp rise of soft gluon density at very small values of x and thus 
would eventually lead to violation of unitary bound on physical cross sections. 

However, the above picture of harder partons spliting into softer ones with 
no further interaction among them is perhaps a bit too naive. Physically we 
know that at very small values of x, the parton density will be high and they 
will spatially overlap. The effects of parton recombination, screening etc., will 
therefore be important. This will eventually lead to saturation of soft gluon 
density at small x. 

It may be more natural to consider soft gluons not as quasi-free partons 
but rather as classical fields with large amplitudes. After summing the effects 
of high density into the clasical fields one can do perturbation theory in the 
background of these large classical fields. 
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2 McLerran-Venugopalan Model 

In the Mclerran-Venugopalan model of small x gluons, one considers a large 
nucleus traveling with very high velocities and thus highly Lorentz contacted. 
It therefore has a high valence quark surface density. This approach can also 
be used for hadrons at sufficiently small values of x where the parton density 
is large. The color charge density p is the only dimentionfuU parameter in the 
problem. Therefore the strong coupling constant a^, as a function of p, will 
be small for large charge densities so that one can use weak coupling meth- 
ods. Since the coupling is weak, the valence quarks do not lose an appreciable 
fraction of their momenta and stay on straight line trajectories and are there- 
fore static sources of color charge. This is the no-recoil approximation. Also, 
due to the high density of glue, one can treat them as classical charges. This 
corresponds to taking a higher dimentional representation of the color algebra. 
Treating the valence charges as classical leads to computing expectation values 
with a Gaussian distribution of the form 

[dp] <^xp{^ j d^xtp^{xt)] (1) 

where /i^ is the average color charge squared per unit area. 

Alternatively, one can solve the Yang-Mills equations in the presence of 
this static color charge. Once a solution is found, one can then compute the 
distribution function with the above Gaussian weight. It was shown inEJ that 
this leads to a distribution function of the Weizsacker- Williams form for soft 
gluons. Quantum corrections to the classical results were also computed and 
similar to standard perturbation theory, it was iound that there are poten- 
tially large logs of ratio of longitudenal momentaa. This led to a longitudcnal 
structure for the color charge density p. Then the problem to consider was to 
solve the classical equation of motion in the presence of this rapidity dependant 
charge density. 

Working in light-cone gauge = and using light cone notation x^ ^ 
(t ± z) and y ^ ln(a;~), we have 

D~A,=g''p{y,xt). (2) 

One can write a formal solution to this equation in terms of p and then 
compute the corelation function ^< A'^Af >p where now both p and 
depend on y. The result is 

= - (-?A|,,)4f'x(..g^).?] (3) 
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where 



(4) 



is the average color charge squared per unit area at rapidity y and resolution 
scale Q^. 



3 RG Equation for x 

Next we need to determine x(y, Q^)- We start with the valence charge density 
and then include hard gluons into the valence charge step by step. In other 
words, we integrate out higher x gluons perturbatively. This generates an effec- 
tive Lagrangian and renormalizes the charge density. Iterating this procedure 
one can derive a RG equation for x- Diagramatically this can be represented 
as 



d X 

dy dQ' 



Figure 1 : diagrams contributing to running of X 

In the first diagram, the propagator of the hard field is in the background 
of classical fields. This propagator was computed incl. In the second diagram 
the hard field propagates in the background of soft modes. This is where 
Sudakov and non-Sudakov form factors are expected to show up. There are 
also diagrams where the hard and soft modes are mixed. 

As a first approxiamation, we will use the following ladder diagram on the 
right hand side of our RG equation where the vertices are eikonalized. It is 
know that ln(.T) terms come from eikonal vertices. 



Figure 2 : eikonalized ladder diagram 

Solid and dashed lines are hard and soft gluons respectively. With this 
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approximation, the RG equation becomes 



dydQ-' 



Gtr{y,Q'). (5) 



It should be emphasized that the right hand side is still a non-linear functional 

ofx(y,Q^). 



4 Solution to RG Equation 

In the high region {Q^ 3> a^x); it can be shown that our RG equation 
reduces to the standard DGLAP evolution equation. Assuming a to be inde- 
pendent of Q^, the approximate solution is 




In the intermediate Q"^ region (Q^ ~ a^x)) the effects of the background fields 
as well as the soft fields arc expected to be important quantitatively. However, 
here we are interested in the qualitative behavior of our RGE. Assuming that, 
as in BFKL, the transverse phase space is a constant, we get the following 
BFKL like behavior 

X ^ ea;p(#— — y). (7) 

Finally in the low region (Q^ <C o?x)^ our RG equation saturates and 
the solution is of the form 

X = Xo + K(yo - y)Q'' (8) 

where k is a slowly varying function of . This leads to saturation of the 
gluon distribution function at small x. As a result, cross sections computed 
with this distribution function would satisfy unitarity bounds. To illustrate 
this saturation, let us consider the behavior of xiv^Q^) some fixed as 
X decreases. In the high transverse momentum region (fcf ^ Q^x)) the dis- 
tribution function will grow as some power of x. As x gets smaller, a^x will 
eventually become larger than and we will be in the saturation region (the 
region between ^^'^ Q;^x) where the distribution function is a slowly 

varying function of k^. This is shown below. 
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Notice that this saturation is due to shrinking of the transverse phase space 
as indicated above. 

5 Acknowledgments 

The work presented here has been done in collaboration with Alex Kovncr, 
Larry McLerran ajid Heribert Weigert. For a full derivation of the results 
presented here seeQ. 

1. L. McLerran and R. Venugopalan,/'/i?/s. i?ew. D49, 2233 (1994); D49, 
3352(1994); D50, 2225(1995). 

2. A. Ayala, J. Jalilian-Marian, L. McLerran and R. Venugopalan, Phys. 
i?eu. D53, 458 (1996). 

3. A. Ayala, J. Jalilian-Marian, L. McLerran and R. Venugopalan, Phys. 
i?et;.D52, 2935 (1995). 

4. J. Jalilian-Marian, A. Kovner, L. McLerran and H. Weigert, hep - 
ph/9606337. 



5 



